Numerical precision radiative corrections to the Dalitz plot of light and 
heavy quark unpolarized baryon semileptonic decays. The cases E+eF 

and A+ Ae+z/ 

Ruben Flores-Mendieta 
Instituto de Fisica, Universidad Autdnoma de San Luis Potosi, 

■ Alvaro Obregon 64, Zona Centra, San Luis Potosi, S.L.P. 78000, Mexico 

o ■ 

O 

■ J. J. Torres 

■ Escuela Superior de Computo del IPN, 



Apartado Postal 75-702, Mexico, D.F. 07738, Mexico 



M. Neri, A. Martinez 



^ I Escuela Superior de Fisica y Matemdticas del IPN, 
O 

^ I Apartado Postal 75-702, Mexico, D.F. 07738, Mexico 

P ■ A. Garcia 



p i' Departamento de Fisica, Centro de Investigacion y de Estudios Avanzados del IPN, 



Apartado Postal 14-740, Mexico, D.F. 07000, Mexico 



(Dated: February 2, 2008) 

[xj ■ Abstract 

We propose and discuss a numerical use of our previous precision results for the radiative corrections to 
unpolarized spin one-half baryon semileptonic decays, which is not compromised to fixing the form factors 
at prescribed values. We present various crosschecks and comparisons with other results available in the 
literature of such analytical radiative corrections. Our analysis, however, is general and applies to all charge 
assignments to the baryons allowed by heavy quarks. The procedure is exemplified with the processes 
rP Y.+eV and A+ Ae+i/. 



PACS numbers: 14.20.Lq, 13.30.Ce, 13.40.Ks 
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I. INTRODUCTION 



The analysis of high statistics experiments of spin one-half baryon semileptonic decays requires 
the inclusion of precision radiative corrections (RC). Despite the important progress achieved in 
the understanding of the fundamental interactions with the standard model iJ], the calculation of 
RC is still committed to model dependence and so are the experimental analyses which use these 
calculations. Fortunately, up to order (a;/7r)(g/Mi) one can show |Q,|3iy,Q] that all the model 
dependence of RC can be absorbed into the already present form factors and that the remaining part 
is model independent. The measurement of such effective form factors is then model independent. 
The results obtained for the RC can be extended to cover all the six possible charge assignments 
of the baryons allowed when heavy quarks are involved |Q]. 

Although one can assert that the RC can be calculated reliably and once and for all to a high 
degree of precision, there still remains the question of their practical use in an actual Monte Carlo 
simulation. It is this question which we shall address in this paper. In practice it is convenient to 
introduce RC into the simulation in a numerical form. However, our previous results are analytical 
l^lsl. Other results in numerical form found in the literature O] are committed to fixing the form 
factors at given values. This last may be not too bad an approximation in hyperon semileptonic 



decays, but it is unacceptable when heavy quarks are involved because the Cabibbo theory |l|] no 
longer provides a reliable guidance. It is therefore desirable to be able to produce numerical RC 
which allow the form factors to be varied freely in the simulation. The goal of this paper is to 
propose a method to do this last and also to illustrate it. 

In Sec. II we briefly review our previous results for the RC of a neutral decaying baryon (NDB) 
and of a negatively charged decaying baryon (CDB). We shall not reproduce the long detailed 
analytical results. They can be found in Refs. yi |5|]. We shall only reproduce in Appendix A 
a convenient rearrangement of several terms which cannot be found in those references. Before 
proceeding towards our goal, in Sec. Ill we shall present three tables with numerical crosschecks 
of the results of Refs. [4J and [5] and a comparison with the results of Ref. ||2|]. The reasons for 
doing this are to help the user gain confidence in the so long analytical results and to provide 
him/her with numbers which serve as reference. In Sec. IV we discuss how to produce arrays of 
numerical RC which keep the form factors uncompromising. As an illustration and also to provide 
the user with reference numbers we shall produce two tables with numerical arrays. Section V is 
reserved for a summary and conclusions. 
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II. RADIATIVE CORRECTIONS TO THE DALITZ PLOT 



First we recall our notation and conventions and next we briefly review the RC to the Dalitz 
plot of unpolarized baryon semileptonic decays. 

The uncorrected transition amplitude for the decay A BlVi is 

G 

Mo = -^^b[P2)W^,{pi,P2)ua{piWi{1)0^v^{p^)], (1) 

where 



+ 



75- (2) 



Here = 7/^(1 + 75), ? = Pi — P2 is the four- momentum transfer, /j and gi are the Dirac vector 
and axial-vector form factors, respectively, and our conventions for the 7-matrices are those of 
Refs. flisl. All the form factors are assumed to be real in this work. In addition to the virtual 
RC to A ^ BlVi one must include the bremsstrahlung RC from the decay A BWi'j. The four- 
momenta and masses of the particles involved in these processes will be denoted by pi = (Ei, pi) 
and Ml, p2 = {E2, P2) and M2, / = {E, 1) and m, pi, = {E^, p„) and rriu = 0, and k = {lu, k) and 
= 0, respectively. The direction of a vector p will be denoted by the unit vector p. p2, 1, and 
p„ will also denote the magnitudes of the corresponding three-momenta when we specialize our 
calculation to the center-of-mass frame of A. No confusion is expected because in this situation 
our expressions will not be manifestly covariant. When a real photon 7 is emitted, the neutrino 
energy becomes E^ = E^ — uj. 

We shall assume that, even if there is no provision to detect real photons, these latter can be 
experimentally discriminated by energy and momentum conservation. Therefore, our expressions 
will be limited to the RC to the three body region of the Dalitz plot of the uncorrected decay 1^0]. 

The complete expression for the Dalitz plot with virtual and bremsstrahlung RC containing 
terms up to order {a/T:){q/ Mi) can be compactly expressed as 



dVi = dfl 



A', + -eu 

71 



(3) 



The phase space factor is d(] = {G^ /2){dE2dEdVL2dipi/ {27:f]2Mi. The index z = G,N 
covers the CDB A- B^l'Vi and the NDB A° B+l'Vi cases, respectively. The full RC 
are contained in the functions Qu. They are given either in terms of triple integrations over the 



photon variables, to be performed numerically, or in an analytical form, where all such integrations 
have been performed explicitly. We shall not reproduce all the long expressions contained in 9j/. 
Instead we shall give the referencing necessary to find their analytical detailed expressions in our 
previous work. We shall make one exception: some expressions can be reduced substantially. For 
the sake of completeness, these expressions are given explicitly in Appendix A. 
Still in a compact form, the RC for the CDB case are 

Gci = B[{<i>c + Ico) + B'^,<^'c + C'a, (4) 

B[, <^'c, and are found in Eqs. (11), (6), (7) and (12) of Ref. fl. Icq is found in Eq. (52) 
of Ref. jsl. 0(7/ is infrared finite, the infrared divergence cutoff of the virtual RC contained in $c 
is cancelled away by its counterpart in the Ico of the bremsstrahlung RC. Here it is important to 
remark that Icq is equivalent to the Io{a) of Eq. (26) plus the first and second summands of Eq. (27) 
in Ref. U. Thus the last summand in Eq. (27) together with Eqs. (28)-(30) of Ref. fl give the 
triple integration form of C'^. Its analytical form is found in Eqs. (43)-(45) of this reference. After 
some algebraic rearrangements, it can be written as 

16 

C'a = HqOq + ^ Hl9i, (5) 

i=2 

the explicit forms of these H'^ are given in Appendix A. H[ does not appear in C^. It is identified 
with the B[ of the first term in Eq. ©, i.e., H[ = B[. 
The compact form of the RC in the NDB case is 

e,v/ = B[{<!>N + Im) + 5^^!$^ + C'^ + C'^A. (6) 

B[ and in this equation are the same of Eq. ©. $jv and are found in Eqs. (7) and (8) of 
Ref. [5], respectively, once the identifications $Ar = 2Re0 and $^ = 2mRe</)' are made. Also in 
this reference the Ino is given in Eq. (40) and B'^^ must be identified with A'^ of Eq. (15). As in 
the charged case, the infrared cutoff contained in $iv and Jtvo has been cancelled away. The triple 

n 

integration form of C'^j^ is found in Eqs. (38), (39), (44), and (45) of Ref. |5]. Its analytical form 
is found in Eqs. (48)-(50) of this reference. However, it can still be reduced into 

C'j^A = DiPni + D2PN2, (7) 
where Di = f'^ + "ig'^ , D2 = f'^ — g'^, and the functions p^i and piV2 are given in Appendix A. 



4 



Expressions Q and Q are model independent, finite in the infrared and ultraviolet limits, and 
are valid to order {a/ii){q/Mi). All the model dependence amounts to six constants which have 
been absorbed in the form factors of Eq. @. This fact is indicated by putting primes on such form 
factors wherever they appear. Thus, A'^ of Eq. @ has the same form as in the uncorrected case. 
Explicitly, it is found in Eq. (10) of Ref. \M- 



ni. CROSSCHECKS 



Before proceeding towards our goal it is convenient to help the reader gain confidence in our 
previous analytical results, as well as to compare them with other results available in the literature. 

Although the triple integration form of the RC in Eq. © is not practical for its use in a Monte 
Carlo simulation, it is very useful to crosscheck the so long and tedious analytical form of such 
corrections. For this latter purpose, we make numerical comparisons of both forms at fixed values 
of E and E2 over the Dalitz plot. The form factors must also be fixed at predetermined values. 
All these crosschecks were satisfactory and it is not necessary to display all the details here. 
Accordingly, we shall present a minimum of illustrative cases. We shall discuss the processes 
S~ neV and H° S+ez7 as examples of CDB and NDB decays. 

The results are displayed in Tables I and II. The form factors, at zero-momentum transfer, have 
been given the arbitrary values /i = 1.0, /2 = —0.97, /a = 0.789, gi = —0.34, g2 = —1.567, and 
^3 = 0.766 for S- ^ neV and /i = 1.0, /a = 1.853, /g = -0.432, gi = 1.267, ^2 = -0.768, 
and (73 = 1.765 in S° — > S+ez7. The anomalous magnetic moments of the baryons contribute 
to the RC and we use k(S-) = 0.3764M^, = -0.6661Afjv, f^in) = 1.9130M^, and 

= O.8895M7V, where Mn is the nuclear magneton. These values were extracted from the 
total magnetic moments given in Ref. using Eq. (22) of Ref. 1I4]. The anomalous magnetic 
moment of the electron is neglected due to the smallness of its contribution. The values of the 
masses also come from Ref. yj]. 

The triple numerical integration results are displayed in entries (a) and the analytical results are 
displayed in entries (b) of those two tables. The energies E and E2 enter through S = E / E^ and 
cr = E2/M1. Em, a"^""^, and o""^™ are determined using the boundaries of the three body region 
given in Ref. iQ]. Throughout these tables one can appreciate a satisfactory agreement within two 
decimal places and the third place being close within rounding effects. 

We have performed another comparison with the numerical results of Ref. lltzh . Again we shall 
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a (a) 



0.8077 


0.1218 


0.1328 0.0963 0.0421 


-0.0159 


-0.0681 


-0.1058 


-0.1212 


-0.1063 


-0.0506 


0.8056 0.1982 


0.2033 0.1443 0.0645 


-0.0166 


-0.0856 


-0.1317 


-0.1452 


-0.1170 


-0.0374 


0.8035 




0.2043 0.1486 0.0698 


-0.0107 


-0.0785 


-0.1223 


-0.1322 


-0.0991 


-0.0141 


0.8014 




0.1974 0.1476 0.0725 


-0.0050 


-0.0699 


-0.1104 


-0.1164 


-0.0791 




0.7993 




0.1871 0.1443 0.0738 


0.0003 


-0.0609 


-0.0976 


-0.0996 


-0.0582 




0.7972 




0.1396 0.0743 


0.0051 


-0.0519 


-0.0844 


-0.0822 


-0.0366 




0.7951 




0.1339 0.0741 


0.0097 


-0.0429 


-0.0710 


-0.0643 


-0.0143 




0.7930 




0.1276 0.0734 


0.0138 


-0.0340 


-0.0573 


-0.0458 






0.7909 




0.0721 


0.0176 


-0.0251 


-0.0432 


-0.0263 






0.7888 




0.0700 


0.0209 


-0.0163 


-0.0284 


-0.0051 






0.7867 






0.0233 


-0.0075 


-0.0121 
















(b) 










0.8077 


0.1218 


0.1328 0.0963 0.0421 


-0.0157 


-0.0672 


-0.1042 


-0.1194 


-0.1049 


-0.0502 


0.8056 0.1982 


0.2033 0.1443 0.0645 


-0.0163 


-0.0847 


-0.1302 


-0.1436 


-0.1159 


-0.0371 


0.8035 




0.2043 0.1486 0.0698 


-0.0105 


-0.0778 


-0.1211 


-0.1308 


-0.0981 


-0.0140 


0.8014 




0.1974 0.1476 0.0725 


-0.0049 


-0.0692 


-0.1092 


-0.1152 


-0.0783 




0.7993 




1 871 1443 0738 


0.0004 


—0.0603 


—0.0966 


—0.0986 


—0.0576 




0.7972 




0.1396 0.0743 


0.0052 


-0.0514 


-0.0836 


-0.0814 


-0.0363 




0.7951 




0.1339 0.0741 


0.0097 


-0.0425 


-0.0703 


-0.0637 


-0.0142 




0.7930 




0.1275 0.0733 


0.0138 


-0.0337 


-0.0568 


-0.0453 






0.7909 




0.0720 


0.0176 


-0.0249 


-0.0428 


-0.0260 






0.7888 




0.0700 


0.0209 


-0.0162 


-0.0282 


-0.0050 






0.7867 






0.0232 


-0.0074 


-0.0120 









(5 0.0500 0.1500 0.2500 0.3500 0.4500 0.5500 0.6500 0.7500 0.8500 0.9500 

^TTiax 0.8078 0.8078 0.8078 0.8078 0.8078 0.8078 0.8078 0.8078 0.8078 0.8078 
o-""" 0.8043 0.7978 0.7925 0.7884 0.7857 0.7847 0.7854 0.7884 0.7939 0.8023 

TABLE I: Values of in S~ neV decay by (a) integrating it numerically and (b) evaluating it 
analytically. is given in units of GeV^. 

limit ourselves to a minimum of examples, since these comparisons were also quite satisfactory. 
For definiteness, we present the case of H~ KeV. For this comparison we must use the values 
of the form factors and of the energies chosen in this reference. Thus /i = 1.0, gi = 0.249, /2 = 
—0.065, and 92 = fs = 93 = 0.0. The results are displayed in Table III. The entries correspond 
to the relative correction defined as Cr = 100{dT — dTQ)/dTQ, where dV is given by Eq. Q 
with i = C and dT^ is the uncorrected decay rate. The upper entries (a) were calculated with our 
analytical expression, using k(E!^) = — 0.5940MAr, and the lower ones (b) come from Ref. 
The agreement is quite reasonable and the small differences one observes can be explained by 
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0.0683 0.0702 0.0371 


-0.0112 


-0.0615 


-0.1001 


-0.1073 


-0.0439 




0.9066 




0.0749 0.0429 


-0.0057 


-0.0553 


-0.0899 


-0.0860 






0.9062 




0.0794 0.0496 


0.0019 


-0.0457 


-0.0744 


-0.0578 






0.9058 




0.0568 


0.0110 


-0.0332 


-0.0547 


-0.0240 






0.9054 




0.0639 


0.0208 


-0.0189 


-0.0316 








0.9050 






0.0293 


-0.0045 


-0.0068 









(5 0.0500 0.1500 0.2500 0.3500 0.4500 0.5500 0.6500 0.7500 0.8500 0.9500 

^TTiax 0.9091 0.9091 0.9091 0.9091 0.9091 0.9091 0.9091 0.9091 0.9091 0.9091 
o-""" 0.9083 0.9070 0.9059 0.9051 0.9047 0.9046 0.9049 0.9055 0.9066 0.9082 

TABLE II: Values of C^^ in rP S+eI7 decay by (a) integrating numerically and (b) evaluating the 
analytic expression. is given in units of GeV^ and is multiplied by 100. 

the differences in assumptions. In Ref. in the bremsstrahlung RC the baryons were assumed 
point-like and higher (a;/7r)(g/Mi)" contributions [n > 2) were included. Instead, we used the 
theorem of Low and kept only the model independent terms. 

Apart from illustration purposes, Tables I-III provide the user of our RC results with numbers 
to compare with. 
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a (a) 



0.8558 


12.4 


4.6 


2.5 


1.1 


0.0 


-1.0 


-2.1 


-3.4 


-5.1 


-8.5 


0.8546 


60.4 


6.2 


3.3 


1.7 


0.4 


-0.7 


-1.9 


-3.3 


-5.2 


-10.0 


0.8534 




8.1 


3.8 


1.9 


0.5 


-0.7 


-2.0 


-3.5 


-5.6 




0.8522 




12.5 


4.4 


2.1 


0.6 


-0.7 


-2.1 


-3.7 


-6.2 




0.8510 




42.9 


5.3 


2.3 


0.6 


-0.8 


-2.3 


-4.0 


-7.1 




0.8498 






7.2 


2.6 


0.6 


-0.9 


-2.5 


-4.4 


-10.1 




0.8485 






13.0 


3.1 


0.7 


-1.0 


-2.8 


-5.1 






0.8473 








4.2 


0.8 


-1.2 


-3.2 


-6.6 






0.8461 








7.6 


0.9 


-1.5 


-4.0 








0.8449 










1.3 


-2.3 


-12.4 


















(b) 












0.8558 


12.4 


4.6 


2.5 


1.1 


0.0 


-1.0 


-2.1 


-3.4 


-5.1 


-8.5 


0.8546 


60.5 


6.2 


3.3 


1.7 


0.4 


-0.8 


-1.9 


-3.3 


-5.3 


-10.0 


0.8534 




8.1 


3.8 


1.9 


0.5 


-0.7 


-2.0 


-3.5 


-5.6 










4 4 


2 1 


6 


—0 8 


—2 1 


—3 7 


— o.z 




0.8510 




42.9 


5.3 


2.3 


0.6 


-0.8 


-2.3 


-4.0 


-7.1 




0.8498 






7.2 


2.6 


0.6 


-0.9 


-2.5 


-4.4 


-10.1 




0.8485 






13.0 


3.1 


0.7 


-1.0 


-2.8 


-5.1 






0.8473 








4.2 


0.8 


-1.2 


-3.2 


-6.6 






0.8461 








7.6 


0.9 


-1.5 


-4.0 








0.8449 










1.3 


-2.3 


-12.5 









5 0.05 0.15 0.25 0.35 0.45 0.55 0.65 0.75 0.85 0.95 

^max 0.8565 0.8565 0.8565 0.8565 0.8565 0.8565 0.8565 0.8565 0.8565 0.8565 
^min 0.8545 0.8510 0.8482 0.8461 0.8448 0.8444 0.8450 0.8466 0.8495 0.8538 

TABLE III: Values of Cr in KeV decay, (a) corresponds to our analytical expression and (b) 

corresponds to Ref. jtzl. 

IV. NUMERICAL FORM OF THE RADIATIVE CORRECTIONS 

As we have seen in the last section, the numerical RC we have discussed are committed to 
fixing the several form factors at prescribed values. This is highly undesirable in an experimental 
analysis. Clearly on the one hand, in the minimization of such analysis, it is necessary to allow the 
form factors to be varied freely and, on the other hand, it is convenient to use the RC in a numerical 
form. In this section we shall discuss a procedure to obtain numerical RC that are not committed 
to fixed values of the form factors and whose use in a simulation amounts to a form analogous to 
matrix multiplication. 

The RC in Eq. Q are quadratic functions of the form factors, so in general they can be 
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expressed as 

6 

Qm=Yl <f^fr (8) 
i<j=i 

In this equation we have momentarily changed our notation and redefined gi = fi,g2 = h, and 
fi'3 = /e; the restriction i < j reduces the sum to 21 terms. The coefficients a™ are functions of E 
and E2 over the Dalitz plot. The index m takes the values m = CI, NI. 

One can calculate the coefficients at fixed points (E, E2) using the analytical results of 
Refs. 141] and and allow such points to cover a lattice over the Dalitz plot. This lattice should 
match at least the bins defined in the experimental arrangement, although for precision RC it 
should probably be made larger so as to allow several {E, E2) points within each bin. 

To calculate the coefficients it is not necessary to rearrange our final results, either analytical 
or to be integrated, so that they take the form dH). One can calculate them following a systematic 
procedure. One chooses fixed {E, E2) points. Then one fixes /i = 1 and fi = 0, i ^ 1 and obtains 
a^; one repeats this calculation for /2 = 1 and /j = 0, i 7^ 2 to obtain a^, and again until /e = 1, 
/j = 0, i 7^ 6 and are obtained. Next, one repeats the calculation with /i = 1, /2 = 1, = 
i 7^ 1, 2 and from this result one subtracts and a^; this way one obtains the coefficient a^. 
One repeats this last step changing i and j until all the interference coefficients a^, i j, have 
been calculated. 

To illustrate all this and to further discuss it we have produced arrays presented in two tables, 
selecting in each one ten points (E, E2) over the Dalitz plot. We have chosen two examples, 
At — > Ae+z/ of a CDB case which is displayed in Table IV and H° S+eI7 of a NDB case 
llOll which is displayed in Table V. The former also serves as an example of a heavy quark decay. 
As in the previous section, the more important purpose is to provide the user with numbers to 
compare with. The arrays of these two tables were obtained using the RC in the analytical form. 
In order to help the interested reader reproduce the entries of Tables IV and V, the (5, o) points 
were chosen such that a = M2/M1 + jAa, with j = 1, 2, . . . , n and Aa = {E^ - M2)/{nMi), 
with = (l/2){(Mi - mY + M|}/(Mi - m). Besides, in the A+ case we used the formulas 
for the charged assignments , B^, of the CDB case of the previous sections and then applied 
the rules of Ref. 1^ to obtain the results for the charge assignments A+, B^, /+ of this particular 
case. In these tables we have restored our standard notation for the axial-vector form factors gi, 
g2, and g^. The masses used are those of Sec. Ill, comes from Ref. 1 1], and we assume an 
estimate for k(A+) = O.IIO6M7V. 
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The first fact that appears in these Tables IV and V is that the RC do not depend on the form 
factor products figs, /25'3, f3gi, 1^92, and f^g^. The nonappearance of these products cannot 
be seen easily in our final results of Sec. IT. The other fact is that the nonzero RC to each form 
factor product vary appreciably from one {E, E2) point to another. This means that replacing the 
precision results of Sec. II with an array of only a few columns over the Dalitz plot is far from 
satisfactory. Therefore the lattice of {E, E2) points must be much finer than only a few points. 

The use of this third presentation of RC is very practical in the sense that such RC can be 
calculated separately and only the arrays should be feed into the Monte Carlo simulation. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have proposed and discussed a numerical form of the model independent RC 
to baryon semileptonic decays which is not compromised to fixing the values of the form factors 
at prescribed values. This last would not be acceptable in decays involving heavy quarks. The 
numerical RC are organized in arrays. Each column of the array is calculated at fixed (£■, £'2) and 
applied as a matrix multiplication to the form factor products. The £'2) points cover a lattice 
over the Dalitz plot. 

The inclusion of precision RC in an experimental analysis is reduced to simple multiplication, 
avoiding the calculation of the long analytical expressions as E and E2 are varied. However, in 
a precision experiment possibly involving 150, 200, and even 300 bins over the Dalitz plot the 
number of columns in the RC arrays should be at least just as many. It may be required that 
several columns be reproduced in finer subdivisions within each bin, possibly four, eight, or even 
more. For example, one may require that the numerical changes of the coefficients between 
neighboring [E, E2) points do not exceed two decimal places within rounding of the third decimal 
place. 

To close, let us recall that our results are general within our approximation. They can be applied 
in the other four charge assignments of baryons involving heavy quarks and whether the charged 
lepton is e^, or t^. They are model independent and not compromised to fixing the form 
factor at determined values. Let us finally emphasize that Tables I-V should serve for comparison 

purposes and help the user to check his/her use of our analytical results. 
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(0.15,0.5995) (0.45,0.5995) (0.75,0.5995) (0.95,0.5995) (0.25,0.5602) (0.55,0.5602) (0.85,0.5602) (0.45,0.5210) (0.75,0.5210) (0.65,0.4948) 



f! 


4.639 




5.764x10" 


-1 


-4.632 




-2.647 




3.369 




-7.737x10" 


-1 


-3.403 


1.176 




-2.016 




-3.654x10" 


-1 


fi 


6.914x10" 


-1 


1.007x10" 


-1 


-6.431x10" 


-1 


-7.604x10" 


-1 


5.006x10" 


-1 


-6.431x10" 


-2 


-7.948x10-1 


8.519x10" 


-2 


-3.917x10- 


1 


-7.668x10" 


-2 


fi 


-2.944x10" 


-1 


-2.499x10" 


'2 


3.269x10" 


-1 


-6.177x10" 


-2 


-1.813x10" 


-1 


7.650x10- 


-2 


2.357x10-2 


-1.037x10" 


-1 


6.046x10- 


■2 


1.093x10" 


-2 


a! 


5.963 




8.010x10- 


-1 


-5.674 




-6.435 




6.416 




-1.008 




-1.013x10+1 


2.271 




-1.038x10+ 


■1 


-9.480 




al 


4.804x10" 


- 1 


8.761x10" 


-2 


-3.905x10" 


- 1 


-9.652x10" 


- 1 


7.038x10- 


- 1 


-1.568x10" 


-2 


-1.565 


2.244x10" 


- 1 


-1.971 




-2.324 




al 


-3.010x10" 


-2 


-2.555x10" 


-3 


3.343x10" 


-2 


-6.316x10" 


-3 


-1.244x10" 


-2 


5.250x10" 


-3 


1.618x10-3 


-3.359x10" 


-3 


1.959x10- 


3 


7.272x10" 


-5 


fifi 


1.011 




1.398x10" 


-1 


-9.471x10" 


-1 


-1.170 




1.055 




-1.736x10" 


-1 


-1.517 


3.411x10" 


-1 


-9.479x10- 


1 


-1.988x10" 


-1 


fifi 


5.303x10" 


-1 


3.216x10" 


-2 


-6.074x10" 


-1 


1.133x10" 


-1 


3.281x10" 


-1 


-1.603x10" 


-1 


-4.672x10-2 


1.912x10" 


-1 


-1.231x10- 


1 


-2.662x10" 


-2 


f2h 


2.176x10- 


-1 


1.947 xlO" 


-2 


-2.403 xlO" 


-1 


4.551x10" 


-2 


1.561x10" 


1 


-6.085x10" 


-2 


-1.937x10^2 


1.038x10" 


-1 


-5.612x10- 


2 


-9.102x10" 


-3 


9192 


-2.084 




-3.016x10- 


-1 


1.792 




4.204 




-3.801 




3.736x10" 


-1 


7.428 


-1.481 




8.939 




9.367 




9193 


-5.523x10" 


-1 


-5.975x10" 


-2 


5.949x10" 


-1 


-1.138x10" 


-1 


-3.635x10" 


-1 


1.315x10" 


-1 


4.321x10-2 


-2.196x10" 


-1 


1.165x10- 


1 


1.785x10" 


-2 


9293 


2.176x10" 


-1 


1.947x10" 


-2 


-2.403x10" 


-1 


4.551x10" 


-2 


1.561x10" 


-1 


-6.085x10" 


-2 


-1.937x10-2 


1.038x10" 


-1 


-5.612x10- 


2 


-9.102x10" 


-3 


/l9l 


-6.338x10" 


-3 


-7.034x10" 


-1 


-6.223x10" 


-2 


2.432 




-8.672x10" 


-1 


-1.921 




2.286 


-2.952 




3.431x10- 


■1 


-1.730 




/l92 


-6.962x10" 


-3 


3.520x10" 


-1 


2.672x10" 


-2 


-1.246 




4.037x10" 


-1 


9.504x10" 


-1 


-1.183 


1.440 




-2.088x10- 


■1 


8.390x10- 


-1 


/l93 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


0.000 




0.000 




0.000 




/291 


-2.744x10" 


-2 


-1.059 




-9.851x10" 


-2 


3.618 




-1.348 




-2.898 




3.387 


-4.477 




4.752x10- 


1 


-2.625 




/292 


-1.036x10" 


-2 


5.239x10" 


-1 


3.977x10" 


-2 


-1.855 




6.009x10" 


-1 


1.415 




-1.761 


2.143 




-3.107x10- 


1 


1.249 




/293 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


0.000 




0.000 




0.000 




/391 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


0.000 




0.000 




0.000 




/392 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


0.000 




0.000 




0.000 




/393 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 


0.000 




0.000 




0.000 





TABLE IV: Numerical arrays of the coefficients afj^ in GeV^ of Eq. (8) evaluated at ten points {S, a) (headings of columns) over the unpolarized Dalitz 
plot of A+ —>■ Ke^v decay. 



(0.05,0.9086) (0.35,0.9086) (0.65,0.9086) (0.95,0.9086) (0.25,0.9073) (0.55,0.9073) (0.75,0.9073) (0.45,0.9059) (0.65,0.9059) (0.55,0.9050) 



ft 


4.874x10" 


-2 


9.845x10" 


-2 


1.783x10" 


-2 


—1.968x10' 


_2 


6.750x10' 


-2 


2.728x10" 


-2 


— 1.885x10" 


-2 


2.339x10" 


-2 


—8.097x10" 


-3 


—1.071x10" 


-3 


/I 


5.047x10" 


-4 


8.128x10" 


-4 


1.708x10" 


-4 


-2.698x10' 


-4 


6.142x10' 


-4 


2.818x10" 


-4 


—9.249x10" 


_5 


2.133x10" 


-4 


2.666x10" 


-6 


1.750x10" 


-5 


/f 


—2.270x10" 


-4 


—6.500x10" 


-4 


—9.241x10" 


.5 


3.358x10' 


_5 


—2.443x10' 


-4 


—6.334x10" 


.5 


1.487x10" 


-4 


—3.114x10" 


_5 


8.899x10" 


.5 


2.953x10" 


-5 


2 

9i 


8.473 X 10" 


_2 


1.201 X 10" 


_1 


2.695 X 10" 


-2 


-5.073x10' 


_2 


1.718x10' 


_1 


8.107x10" 


-2 


— 1.978x10" 


_2 


1.378x10" 


_1 


7.920x10" 


.3 


3.924x10" 


-2 


2 


3.134x10" 


-4 


1.856x10" 


-4 


8.773x10" 


.5 


— 2.664X 10' 


-4 


7.566x10' 


-4 


4.189X 10" 


-4 


5.428x10" 


_5 


9.168x10" 


-4 


1.964x10" 


-4 


3.790x10" 


-4 


2 


—5.090x10" 


_7 


— 1.457x10" 


-6 


—2.072x10" 


.7 


7.528x10' 


-8 


—3.654x10' 


_7 


—9.475x10" 


-8 


2.225x10" 


_7 


—2.331x10" 


-8 


6.661x10" 


-8 


7.372x10" 


-9 


/1/2 


5.638x10" 


-4 


9.122x10" 


-4 


1.875x10" 


-4 


—3.016x10' 


-4 


8.794x10' 


-4 


3.890x10" 


-4 


— 1.641x10" 


-4 


3.700x10" 


-4 


—4.900x10" 


-5 


8.182x10" 


-6 


/1/3 


2.474x10" 


-4 


7.098x10" 


-4 


9.664x10" 




— 3.723 X 10 


-5 


2.585x10' 


-4 


6.115x10" 


-5 


— 1.671x10" 


-4 


2.655x10" 


-5 


-1.010x10" 


-4 


— 3.519X 10" 


-5 


/2/3 


2.291 X 10 


-5 


6.557 X 10" 


-5 


9.365 X 10 


-5 


— 3.3s2x 10 




2.544X 10 




D.s20x 10 




— 1.51 1 X 10" 


-5 


3.760x10 


-5 


— 9.076x10 


-5 


— 2.900X 10 


-6 


9192 


—7.212x10" 


-3 


—4.911x10" 


-3 


— 1.891x10" 


-3 


5.952X 10' 


-3 


— 2.032X 10' 


_2 


— 1.045x10" 


.2 


2.576x10" 


-4 


—2.196x10" 


-2 


—3.032x10" 


-3 


— 7.708 X 10" 


-3 


9193 


—2.541 X 10 


-4 


— 7.244 X 10 


-4 


— 1.072X 10 


-4 


3.692X10 


-5 


Oil ^iin- 

—2.811 XlO 


-4 


— 7.867X10 


-5 


1.614X 10 


-4 


— 4.222 X 10 


-6 


9.575X10 


-6 


3.032X10 


-5 


9293 


2.291x10" 


-6 


6.557x10" 


-6 


9.365x10" 


-6 


-3.382x10' 


-6 


2.544x10' 


-5 


6.820x10" 


-6 


-1.511x10" 


-6 


3.760x10" 


-6 


-9.076x10" 


-6 


-2.900x10" 


-6 


/191 


-2.678x10" 


-3 


1.609x10" 


-4 


9.352x10" 


.4 


-2.700x10' 


-3 


-2.051x10' 


-3 


3.242x10" 


-3 


1.489x10" 


-3 


3.140x10" 


-3 


2.945x10" 


-3 


3.124x10" 


-3 


/l92 


2.484x10" 


-4 


-1.926x10" 


-5 


-8.964x10" 


-5 


2.580x10' 


-4 


1.858x10' 


-4 


-3.063x10" 


-4 


-1.364x10" 


-4 


-2.898x10" 


-4 


-2.690x10" 


-4 


-2.829x10" 


-4 


/l93 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




/251 


-4.874x10" 


-3 


4.487x10" 


-4 


1.824x10" 


-3 


-5.142x10' 


-3 


-3.363x10' 


-3 


6.360x10- 


-3 


2.878x10" 


-3 


6.339x10" 


-3 


5.704x10" 


-3 


6.106x10- 


-3 


/292 


4.732x10" 


-4 


-3.667x10" 


-6 


-1.707x10" 


-4 


4.914x10' 


-4 


3.539x10' 


-4 


-5.834x10" 


-4 


-2.598x10" 


-4 


-5.520x10" 


-4 


-5.124x10" 


-4 


-5.388x10" 


-4 


/293 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




/391 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




/392 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




/393 


0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 




0.000 





TABLE V: Numerical arrays of the coefficients afj^ in GeV^ of Eq. (8) evaluated at ten points [6, a) (headings of columns) over the unpolarized Dalitz 
plot of S° —>■ S+ez^ decay. 
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APPENDIX A: COLLECTION OF THE m AND p, FUNCTIONS 



In this section we present the H'- and the pi functions of C'^ in Eq. ^ and C'j^^ in Eq. Q, 
respectively. They read 



K=P2l 



Qs + QaEI P2IE 

2 Ml 



2 



^, = P^^C^ f^') [-EliQ^ + Q3) + Q,{E + Elf + {Q2 + Q,)p2lyo + 



E'^.-E 



{E + El+p2(3yo)Qz 



P2lyo+pl-E{E + E^ 



+ Q4 



ElPp2yo - EiE + El 



+ 



"""^^-^^l [(E + El)B' - EB^] + 2El[A- + ,i(/, - ,,)] 



- AEg^{f2 + 92) + 2mJ^{2E + E°)| 



^{QiE + Q,{2El + 3E) - iQ2 + Q4)E[Eil + /3') + E^ 



+ 2Q4El{El-E)-p2lyo]} 



P2IE 
2Mi 



[~{E + El)B~ + EB++p2PyoC-] 



M,^[El + 2E{1 - P')] + M,^/3'E - Elifj2 + 9i92)il + /3xo) 



H' 



+ fracp2lEMi 

+ E%{f2 + 392) + E[{1- P'){fj2 + ^9192 - 9if2) - /1/2 + 9im - 292)]] 
P2I' 



jgi + 3Q3 -Q2{E + El) - Q,{3E + 7El) + A [(£;0 _ j^^j^- ^ ^ ^0^-] 



+ P2I' 



h- 



E. 



2E 



{2E + El) 2El + ^(2/1/2 + 3^71^72 - ^71/2) + Tr9iU2 + 92) 



Ml 



Ml' 



13 



K = [Qi + Q3- iQ2 + Q4){E + , 



^ 4 



E 



(Qi + Q3) ^ + {Q2 + b2/3|/o - + - E{1 - p^)] + g^^ 



E 



P2I 



+ 



2Mi 



4Mi 



M^ — {El + + p2/3yo) - - /32:o)yi(^2 - /2) - £^(1 - 

e 



+ 



2Mi 



H',-^-f [-Qi - ^3 + {Q2 + g4)(i^ + El)] , 



10 



^ (g2 + 5g4) + |^(2B- + c-) 



+ 



P2^^ 

Ml 



_2Mi— - 3Mi — + 3/1/2 + 4(71(72 - 3(7i/2 
e e 



It turns out that = in this rearrangement. However, the indexing of Eq. (5) is made easier 
by keeping this H[-^ expUcitly. 



and 



with 



2;2 



H' 



2/2 



2/2 



Pil 



12 



2 2Mi Ml 



Ml— + ^1(52-/2) 



2/2 



2/2 



13 



■Qi- 



2;2 



2Mi 



-c- - 



pV_ 

Ml 



Ml— + /2(/l+^l) 



14 4 ^ 4Mi^ ^' 



P2I 



P2I 



P2I 



15 



— (g2 + Qa) + -^B- + 

8 4Mi 2Mi 



Ml— + /2(/l-^l) 



16 



P2/3 

4Mi 



-Ml— + ^1(^2-/2) 
e 



/1/2 - 9'ifl'2 ± 2g'i/2, 

(/l±(?lf + 2/i(/3-/2), 
/l -^l+2/l(/3-/2), 

-5fi(«;i + K2) ± /i^i(«;2 - Ki). 



14 



Here, e is the electron charge, /3 = l/E, yo = {{E^)'^ — P — p^} / {2lp2), xq = —{p2yo + l)/E^ and 
Ki and K2 are the anomalous magnetic moments of the decaying and emitted baryons, respectively. 
The Pi functions of Eq. ^ are 



Pm = P 



2Mi 



13 



and 



with 



P2I 
2Mi 



(i5;0 _ E)9o + m\93 - 9,) - 1{E + 2El)9, 



Pp2yo + E + EI + — 



9_j_ 
2 



+ {El -2E)%- 3l%o + PP2{2E - El)9,2 - Wu - % + 



4:E 



The Qi and 9i functions are found in Refs. 
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